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Abstract 



We obtain a representation of Feller's branching diffusion with logistic growth in terms 
of the local times of a reflected Brownian motion H with a drift that is affine linear in 
the local time accumulated by H at its current level. As in the classical Ray-Knight 
representation, the excursions of H are the exploration paths of the trees of descendants 
^ of the ancestors at time t = 0, and the local time of H at height t measures the population 

I— I size at time t (see e.g. [H]). We cope with the dependence in the reproduction by 

introducing a pecking order of individuals: an individual explored at time s and living at 
^ time t = Hs is prone to be killed by any of its contemporaneans that have been explored 

so far. The proof of our main result relies on approximating H with a sequence of Harris 
paths which figure in a Ray-Knight representation of the total mass of a branching 
particle system. We obtain a suitable joint convergence of together with its local 
" I times and with the Girsanov densities that introduce the dependence in the reproduction. 

in 
o 

1 Introduction 

Feller's branching diffusion with logistic growth is is governed by the SDE 
^ dZt = a^/ZtdWt + {eZt - jZ^) dt, Zq = x > 0, (1.1) 

with positive constants a, 9 and 7. It has been studied in detail by Lambert [10], and models 



the evolution of the size of a large population with competition. The diffusion term in (1.1) 
incorporates the individual offspring variance, and the drift term includes a super-criticality 
in the branching that is counteracted by a killing with a rate proportional to the "number of 
pairs of individuals" . 



For = 7 = 0, equation (1.1) is the SDE of Feller's critical branching diffusion with 
variance parameter . In this case, a celebrated theorem due to Ray and Knight says that Z 
has a representation in terms of the local times of reflected Brownian motion. To be specific, 
let H = {Hs)s>o be a Brownian motion on R_|_ with variance parameter reflected at the 
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origin, and for s,t >0 let Ls{t, H) be the (semi-martingale) local time accumulated by H at 
level t up to time s. Define 

:= inf{s > : (cjV4)L,(0, H) > x}. (1.2) 



Then {a /4)Ls^{t, H), t > 0, is a weak solution of (1.1) with = 7 = 0, and is called the 
Ray-Knight representation of Feller's critical branching diffusion. 

The Ray-Knight representation has a beautiful interpretation in an individual-based pic- 
ture. Reflected Brownian motion H = {Hs)s>q arises as a concatenation of excursions, and 
each of these excursions codes a continuum random tree, the genealogical tree of the progeny 
of an individual that was present at time t = 0. The size of this progeny at time t > is 
times the (total) local time spent by this excursion at level t. Starting with mass x at time 
t = amounts to collecting a local time (4/(T^)x of -ff at level 0. The local time of H at level t 
then arises as a sum over the local time of the excursions, just as the state at time t of Feller's 
branching diffusion, Zt, arises as a sum of the masses of countably many families, each of 
which belongs to the progeny of one ancestor that lived at time t = 0. The path {Hs)o<s<Sx 
can be viewed as the exploration path of the genealogical forest arising from the ancestral 
mass X. We will briefly illustrate this in Section [3] along a discrete mass - continuous time 
approximation. For a more detailed explanation and some historical background we refer to 
the survey |15j . 



The motivation of the present paper was the question whether a similar picture is true 



also for (1.1) with strictly positive 9 and 7, and whether also in this case a Ray-Knight 



representation is available for a suitably re-defined dynamics of an exploration process H. At 



first sight this seems prohibiting since the nonlinear term in (1.1) destroys the independence in 
the reproduction. However, it turns out that introducing an order among the individuals helps 
to overcome this hurdle. We will think of the individuals as being arranged "from left to right" , 
and decree that the pairwise fights are always won by the individual "to the left" , and lethal 
for the individual "to the right" . In this way we arrive at a population dynamics which leaves 



the evolution (1.1) of the total mass unchanged, see again the explanation in Section [3j The 
death rate coming from the pairwise fights leads in the exploration process of the genealogical 
forest to a downward drift which is proportional to Ls{Hs, H), that is, proportional to the 
amount of mass seen to the left of the individual encountered at exploration time s (and 
living at real time Hs) - more rigorously, Lr{t,H) denotes the local time accumulated by 
the semimartingale H up to time r at level t, see the beginning of Section [2] for a precise 
definition. As a consequence, those excursions of H which come later in the exploration tend 
to be smaller - the trees to the right are "under attack from those to the left" . 
In quantitative terms, we will consider the stochastic differential equation 

2 1 26 

Hs = -Bs + -Ls{0,H) + ^s-j LriHr,H)dr, s>0, (1.3) 

where B \s & standard Brownian motion. The last two terms are the above described com- 
ponents of the drift in the exploration process, and the term Ls{0, H)/2 takes care of the 
reflection of H at the origin. We will show 



Proposition 1.1. The SDE (1.3) has a unique weak solution. 



Our main result is the 
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Theorem 1. Assume that H solves the SDE (1.3), and let, for x > 0, Sx he defined as in 
(O). Then {a'^ /4)LsAt, H), t > 0, solves (O). 



We will prove Proposition |1.1| by a Girsanov argument, and Theorem [T] along a discrete 
mass-continuous time approximation that is presented in Section [3j In section |4] we take the 
limit in the total population process along the discrete mass approximation. An important 
step in the proof of Theorem [l| and interesting in its own right, is Theorem [2] in Section 
[5j in which we obtain a convergence in distribution, in the case = 7 = 0, of processes 
that approximate reflected Brownian motion, together with their local times (considered as 
random fields indexed by their two parameters, time and level). A similar convergence result 
was proved in [T7] for piecewise linear interpolation of discrete time random walks and their 
local time. The proof of Theorem[T]is completed in Section[6} using again Girsanov's theorem, 
this time for Poisson point processes. In that section we will also prove Theorem [3| which says 
that in the case ^, 7 > 0, the exploration process together with its local time (now evaluated at 
a certain random time, while the parameter for the various levels is varying) converges along 
the discrete mass approximation. Finally, for the convenience of the reader, we collect in an 
Appendix, at the end of this paper, several results from the literature, concerning tightness 
and weak convergence in the space D, the Doleans exponential and "goodness", and the two 
versions of Girsanov's theorem which we need : the one for Brownian motion and the one for 
Poisson point processes. 

When this work was already completed, our attention was drawn by Jean-Frangois Le Gall 
on the article |14j by J. Norris, L.C.G. Rogers and D. Willams, who proved a Ray-Knight 
theorem for a Brownian motion with a "local time drift" , using tools from stochastic analysis, 
in particular the "excursion filtration" . With a similar methodology we were recently able to 
establish another, shorter but less intuitive, proof of the main result of this paper |16] . 



2 Proof of Proposition 1.1 



For abbreviation we write Ls{t) := Ls{t,H) for the local time at level t > 0, accumulated 
by the continuous semi-martingale H up to time s. For the reader's convenience, we recall a 
possible definition (borrowed from [T9], Theorem VI. 1.2) of that quantity : 



L,{t) ■.= 2{Hs-ty 



'^{Hr>t}dHr. 



Following the Girsanov route, let i? be a standard Brownian motion defined on a proba- 
bility space (fi, J^, P) and let H obey 



Hs = -Bs + Ils{0), 
a 2 



(2.1) 



which can be viewed as Skorohod's equation for refiected Brownian motion, see section 3.6.C 

in [8]. Thus the process H is under the probability measure P a reflected Brownian motion 

^2 



1.1 



with variance parameter 4/a . In order to prove the existence part of Proposition 
suffices to construct, by a suitable reweighting of P, a probability measure P under which 



it 



(77 



a 



dr. 



s > 0, 



(2.2) 
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is a standard Brownian motion. Indeed, (2.2) and (2.1) together imply that under P the 
process H solves (1.3) with B instead of B. 

The Girsanov condition which makes the reweighting possible and whose validity we will 
check is 

Eexp [Ms - J(Af) J =1, s > 0, (2.3) 



with M, := /; [f - ^LriHr)] dBr. 



Remark 2.1. If (|2.3|) holds, then there exists a measure P such that for all s > 

1 



dF 

dF 



exp(M, --(M), 



wher e J^s is the a-field generated by {Hr)o<r<s- By Girsanov's theorem (see Proposition 
7.12 in the Appendix below), B defined by (|2.2[) then is a standard Brownian motion under 



P. Thus, under P, H solves (1.3) with B instead of B. This gives the existence part of 
Proposition |1.1[ as soon as we have (2.3). 



A sufficient condition for (2.3) is provided by the following lemma, which is Theorem 1.1, 
chapter 7 (page 152) in [6]. 

Lemma 2.2. Assume that the quadratic variation of the continuous local martingale M is of 
the form {M)s = Jq Rr dr, and that for all s > there exist constants a > and c < oo such 
that 

Eexp(aii,) < c, < r < s. (2.4) 



Then (2.3) is satisfied. 
In our situation, R^ 



cr7 



Lr{Hr)\ , for which (2.4) is implied by the following 



Lemma 2.3. Let H be a Brownian motion on M_|_ reflected at the origin, with variance 
parameter v"^ . Then for all s > there exists a = a{s, v) > and a constant c < oo such that 

E {exp{aLr{Hrf)) < c, < r < s. 

Proof: Together with a simple scaling argument and a desintegration with respect to Hr, this 
is immediate from the following 

Lemma 2.4. Let (3 be a standard Brownian motion starting at 0, and denote by Li{t) the 
local time accumulated by \P\ at position t up to time 1. There exist constants a > and 
c> ( not depending on t) such that 



E[e''^iW'| = i] < c, i>0. 



(2.5) 



Proof: Denote by Ki{x) the local time of (3 accumulated up to time 1 at position x. First 
observe that for t > 

Li{t) = Ki{t) + Ki{-t) a.s. (2.6) 



For deriving (2.5), by symmetry it suffices to condition under the event = t}. Writing P^ 

uali 

1/2 



for P[. |/3i = x] we conclude from (2.6) and the Cauchy-Schwarz inequality that for all a > 

1/2 



E«[e"^iW ] < (E* 



E' 



(2.7) 
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For u < 1, the distribution of Ki(t) under P* and conditioned under the event that /3 hits t 
first at time u, equals the distribution of \/l — uKi{Q) under F^. Similarly, for ui,U2 < 1, 
the distribution of Ki{—t) under P* and conditioned under the event that {f3v)o<v<i hits —t 
first at time ui and last at time U2, equals the the distribution of y/u2 — ui i^i(O) under F^. 
Consequently, 

\4aXi(0)2l 



E* 



< E^ 



(2i 



a Raleigh distribution. 



By a result due to Levy (see formula (11) in [18j), Ki{0) has under 
i.e. 

P°(i^i(0) >£) = e-i^'. 

This means that Kf{0) is exponentially distributed, and hence, for suitably small S > 0, 



EO 



is finite. Now (2.5) follows from (2.7) and (2.8). 



As stated above. Lemmas 2.2 and 2.3 give together with Remark 2.1 the existence part of 
Proposition 1 1.1[ For its uniqueness part, assume that H is a weak solution of (1.3), governed 
by some measure P. For n S N, we define 



T„ := inf{r > : Lr(i?r) > n}. 



(2.9) 



By a Girsanov transformation we can change the measure P into a measure P under which, 
for aU n G N, 

j-sAnAT, 

BsAnATn + 



AnAT„ 



a 



is a standard Brownian motion stopped at n A T„ 
satisfies 

dHs = -dBs + ldLs{0), 
a I 



dr, s > 0, 
Then, under P, for all n, the process H 

< s < n AT„, 







0. 



Because the weak solution of (2.10) is unique, the law of {HsAnATn)s>o under P is uniquely 
determined, and so is the law of {HsAnAT„)s>o under P, since it can be recovered by the inverse 
Girsanov transformation (note that P and P are mutually absolutely continuous on J-nATn)- 
From Theorem VL1.7, page 225 of [19j one can infer that the local time L = Ls{t) of H is 
jointly continuous on [0, oo) x [0, oo) in its two arguments P-a.s. Consequently, T„ — ?• oo as 
n — )• oo, P-a.s, which shows weak uniqueness for (1.3) and completes the proof of Proposition 

o 



3 A discrete mass approximation 

The aim of this section is to set up a "discrete mass - continuous time" approximation of 
(1.1) and (1.3). This will explain the intuition behind Theorem [T| and also will prepare for 



its proof. 

For X > and G N the approximation of ( 1.1 ) will be given by the total mass Z^'^ of a 



population of individuals, each of which has mass 1/N . The initial mass is Z^'^ = [AxJ/A, 
and Z^'^ follows a Markovian jump dynamics: from its current state k/N , 

yN^.. , {{k + l)/N aiTatekNa^/2 + ke 

Z ' lumps to < „ (3.1) 

[ (A: - 1)/A at rate fcA^o-V2 + A:(A; - l)7/iV. 
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Figure 1: A realization of (the first two trees of) F and (the first two excursions of) its 
exploration path . The i-axis is real time as well as exploration height, the s-axis is 
exploration time. 

For 7 = 0, this is (up to the mass factor l/N) as a Galton- Watson process in continuous 
time: each individual independently spawns a child at rate Na'^ /2 + 9, and dies (childless) at 
rate Na'^/2. For 7 7^ 0, the additional quadratic death rate destroys the independence, and 
hence also the branching property. However, when viewing the individuals alive at time t as 
being arranged "from left to right" , and by decreeing that each of the pairwise fights (which 
happen at rate 27 and always end lethal for one of the two involved individuals) is won by 
the individual to the left, we arrive at the additional death rate 2jCi{t)/N for individual i, 
where Ci{t) denotes the number of individuals living at time t to the left of individual i. 

The just described reproduction dynamics gives rise to a forest F^'^ of trees of descent, 
drawn into the plane as sketched in Figure 1. At any branch point, we imagine the "new" 
branch being placed to the right of the mother branch. Because of the asymmetric killing, 
the trees further to the right have a tendency to stay smaller: they are "under attack" by the 
trees to their left. Note also that, with the above described construction, the F^'^, x > 0, 
are coupled: when x is increased by a new tree is added to the right. We denote the 

union of the F^'^, x > 0, by F^. 

From F we read off a continuous and piecewise linear M+-valued path = (i/f ) 
(called the exploration path of F^) in the following way: 

Starting from the root of the leftmost tree, one goes upwards at speed 2N until one hits 
the top of the first mother branch (this is the leaf marked with D in Figure 1). There one 
turns and goes downwards, again at speed 2N, until arriving at the next branch point (which 
is 53 in Figure 1). From there one goes upwards into the (yet unexplored) next branch, and 
proceeds in a similar fashion until being back at height 0, which means that the exploration 
of the leftmost tree is completed. Then explore the next tree, and so on. For j; > we denote 
by the time at which the exploration of the forest F^'^' is completed. Obviously, for each 
t > 0, the number of branches in F^'^ that are alive at time t equals half the number of 
i-crossings of the exploration path of F^ stopped at . Recalling that the slope of is 
±2A^, we define 

A^ (t) := ^^7^ of i-crossings of between exploration times and s, (3-2) 
where we count a local minimum of at t as two t-crossings, and a local maximum as none. 
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Note that by our convention both s i— )• (t) and t i— )• (t) are right continuous, and in 
particular Aq^ (0) = 0. We cah A^(i) the (unsealed) loeal time of accumulated at height 
t up to time s. This name is justified also by the following oceupation times formula, valid 
for all measurable / : M+ — )• M+ 



f oo 

f{H^) dr = I fit) Af (t) dt, s > 0. (3.3) 
Jo 

The exploration time 5^ which it takes to traverse all of the [Nx\ trees in F^'^ can be 
expressed as 

= mf{s:A^{0)>[Nx\/N}. (3.4) 
Proposition 3.1. The exploration path s i— t- obeys the following stochastic dynamics: 

• At time s = 0, starts at height and with slope 2N. 

• While moves upwards, its slope jumps from 2N to —2N at rate N'^a'^ +A'yN£, where 
£ = (H^) is the local time accumulated by at the current height up to the 
current exploration time s. 

• While H moves downwards, its slope jumps from -2N to 2N at rate N^a"^ + 2Ne. 

• Whenever reaches height 0, it is reflected above 0. 

Proof: We give here an informal proof which contains the essential ideas. (A more formal 
proof can be carried out along the arguments of the proof of Theorem 2.4 in [Ij.) 

Recall that the death rate of an individual i living at real time t is Na"^ /2 + 2jCi{t)/N, 
where Ci{t) is the number of individuals living at time t to the left of individual i. Assume 
the individual i living at time t is explored first at time s, hence Hg = t, and H has slope 2N 
at time s. Because of (3.2), while goes upward, we have Ci{t) = NAs{H^). The rate in t 
is the rate in s multiplied by the factor 2A^ which is the absolute value of the slope. This gives 
the claimed jump rate 2N{Na'^ /2 + 2jAs{H^)) from slope 2A^ to slope —2N, which can be 
seen as the rate at which the "death clock" rings (and leads to a downward jump of the slope) 
along the rising pieces of the exploration path H'^ . On the other hand, the "birth clock" 
rings along the falling pieces of , its rate being Na"^ /2 + in real time and 2N{Na'^ /2 + 6) 
in exploration time, as claimed in the proposition. Note that the process of birth times along 
an individual's lifetime is a homogeneous Poisson process which (in distribution) can as well 
be run backwards from the individual's death time. Also note that, due to the "depth-first- 
search" -construction of , along falling pieces of always yet unexplored parts of the 
forest are visited as far as the birth points are concerned. ■ 

The next statement is a discrete version of Theorem [T| and will later be used for the proof 
of Theorem [T] by taking N ^ oo. 

Corollary 3.2. Let be a stochastic process following the dynamics specified in Proposition 



3.1. and A^ be its local time as defined by (3.2). For x > 0, let be the stopping time 
defined by (3.4). Then 1 1— )■ A^jv(t) follows the jump dynamics (3.1). 



Proof: By Proposition 3.1 , is equal in distribution to the exploration path of the random 
forest F^ . Hence AgN{t) is equal in distribution to Z^'^ , where NZ^'^ is the number of 



branches alive in F^'^ at time t. Since Z^'^ follows the dynamics (|3.1l), so does A 
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The next lemma will also be important in the proof of Theorem [T] 
Lemma 3.3. Let and S!^ be as in Corollary 



3.2 



Then Si 



N 



oo a.s. as x 



oo. 



Proof: Consider x = a/N for a = 1, 2 . . .. Applying (3.3) with s = S^ and / = 1 we obtain 
the equality S"^ = j^K^f^{t)dt. 



3.2 



gN 



(t) follows the jump dynamics (3.1 ), with initial condition 



According to Corollary 
AgN{t) = a/N. By coupling A^^v 
jumps from k/N to {k — 1)/N at rate k{k - 

is stochastically bounded from below by J^^ K^dt, where T2 is the first time at which 
comes down to 2/A^. The latter integral equals a sum of independent exponentially distributed 
random variables with parameters (j — 1)(A^^(T^/2 + 7)), j = 2, . . . ,a. This sum diverges as 



^ (t) with a "pure death process" that starts in a/N and 
l)(iVcrV2 + j/N), we see that A^rf{t)dt 



00. 



4 Convergence of the mass processes Z^'^ as — > 00 

The process {Z^'^, t > 0} with dynamics (3.1) is a Markov process with values in the set 
E]\[ := {k/N, k > 1}, starting from [Nx\/N, with generator A'^ given by 



A^'fiz) = Nz[N"- + 



N 



(4.1) 



+ Nz{N^+^(z-^ 



for any / : E^ — t- M, z G E'tv- (Note that the distinction between symmetric and ordered 
killing is irrelevant here.) Applying successively the above formula to the cases f{z) = z and 
f{z) = z^, we get that 



ZN,X ryN,X . 

t — ^0 + 



N 



z 



■N.x 



z, 



N,x 



f 

Jo 



+ 2 Z 



N,x 



9Z 



N,x 



■N,X I ryN,X 



dr + Ml 

1 



(1) 



N 



a Zj.' + —Z^ + 



N,x I 7 / 



N 



N-^ 



dr 



dr + Ml'^\ 



(4.2) 



(4.3) 



where {M^^\ t > 0} and {M^'-^\ i > 0} are local martingales. It follows from (4.2) and (4.3) 
that 



■t r 



{M(% = [ 
Jo 

We now prove 

Lemma 4.1. For any T > 0, 



2 ryN,X I " yN.X I 

' N 



■N,X I 7 / rp-N 



N 



N 



7N.X 



dr. 



(4.4) 



sup sup E 

7V>10<t<T 



z 



■N.x 



< 00. 
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An immediate Corollary of this Lemma is that {Alj:^^} and {M^'^^} are in fact martingales. 
Proof: The same computation as above, but now with f{z) = z^, gives 



N,x 



N,x 



+ / ^j,(Z^'-)dr + M!^'^ 



t ' 



where {M^^\ t > 0} is a local martingale and for some c > independent of N, 



^iv(^) <c(l + /). 



(4.5) 



(4.6) 



We note that NZ^'^ is bounded by the value at time t of a Yule process (which arises when 
suppressing the deaths), which is a finite sum of mutually independent geometric random 
variables, hence has finite moments of any order. Hence M^'^^ is in fact a martingale. We 
then can take the expectation in (4.5), and deduce from (4.6) and Gronwall's Lemma that 
for < i < r, 



E 

which implies the result. 

We shall also need below the 

Lemma 4.2. For any T > 0, 



7N,X 



< 



z 



N,x 



+ CT 



sup E 

N>1 



sup [Z^ 

0<t<T 



N,x 



< 00. 



Proof: Since from ([42]), zf '"^ < Z^'^ + 9 /g Z^'^'dr + Mi^\ 

sup |Zf < 3|Zo^'"|2 + 3td^ f \Z^'-fdr + 3sup [M^ 

r<t Jo r<t 



1)|2 



This together with (4.4), Doob's L^-inequality for martingales and Lemma 
result. 



4.1 



implies the 



Remark |7.8| in the Appendix combined with (4.2), (4.4) and Lemma 4.1 guarantees that 
the tightness of {Zq^}„>i implies that of {Z^}jv>i in D([0, +00)). 

Standard arguments exploiting (4.2) and (4.3) now allow us to deduce the convergence of 
the mass processes (for a detailed proof, see e.g. Theorem 5.3 p. 23 in |13j). 

Proposition 4.3. As N ^ cc, Z^'^ Z^ , where Z^ is the unique solution of the SDE (1.1) 
and thus is a Feller diffusion with logistic growth. 



5 Convergence of the exploration path in the case ^ = 7 = 0. 



3.1 



with = 7 = 0. The aim of this section 
of H^ and its local time which is 



00 



Let H^ be a stochastic process as in Proposition 
is to provide a version of the joint convergence (as 
suitable for the change of measure that will be carried through in Section [6j This is achieved 
in Theorem [2] and its Corollary 5.1 The proof of Theorem 2 is carried out in two major parts. 
The first part (Proposition 5.3) provides a refined version of the joint convergence of H^ and 
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its local time at level 0, the second part (starting from Lemma 5.5) extends this to the other 
levels as well. 

We define the (scaled) local time accumulated by H^ at level t up to time s as 



Note that this process is neither right- nor left -continuous as a function of s. However since 
the jumps are of size 0{1/N), the limit of as A'^ — )> 00 will turn out to be continuous. In 
fact, we will show that converges as — )• 00 to the semi-martingale local time of the 
limiting process H, hence the scaling factor 4/(T^. It is readily checked that -L^(O) = 
and 



4 



^A^{t), Vs>0,i>0, 



(5.1) 



where A was defined in (3.2). Then with defined in (3.4), we may rewrite 



S!^ = inf{s > : (0) > -^[Nx\/N}. 



(5.2) 



From (5.1) and Corollary 3.2 we see that 



ryN,X 



follows the jump dynamics (3.1 ) in the case ^ = 7 = 0. 

Let {F/, s > 0} be the cadlag {-1, l}-valued process which is such that for a. a. s > 0, 



dH 



N 



ds 



2NV^ . 



We can express in terms of H^^ and V^^ as 



N 



rN 



4 1 
~^2N 



0<r<s 



1 + 



+ 1 



{H^=t}^{Vl 



-1} 



where we put Vq- = +1. Note that any r < s at which f is a local minimum of H counts 
twice in the sum of the last line, while any r < s at which t is a local maximum of H^ is not 
counted in the sum. 
We have 



H. 



N 



N 



2N 



1 + 2 



Vi'dr , 



N 



Na 



(5.3) 



+ — (Lf(0)-L-+(0)) 



where {P^,s > 0} is a Poisson process with intensity N'^a'^. Note that = P^ — N'^a'^s, 
s > 0, is a martingale. The second equation is prescribed by the statement of Proposition 



3.1 (in fact its simplified version in case = 7 = 0) : the initial velocity is positive, whenever 
= 1, it jumps to -1 (i. e. makes a jump of size -2) at rate N'^a'^, whenever = —1, it 
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jumps to +1 (i. e. makes a jump of size +2) at rate N a , and in addition it makes a jump 
of size 2 whenever H'^ hits 0. Writing in the first line of (5.3) as 

l{y^J^=+i} - l{v;^=-i} 

and denoting by the martingale P^ — N'^a'^s, we deduce from (5.3) 

1 



H^ + 



N 



M. 



1,N 



M. 



2,7V 



+ -Lf(0)--L- (0), 



1 



-N 



where 



M. 



1,7V 



N 



and M. 



2,7V 



Na^ Jo 



7V 



(5.4) 



(5.5) 



are two mutually orthogonal martingales. Thanks to an averaging property of the V (see 
step 2 in the proof of Proposition 5.3 below) these two martingales will converge as — t- 00 
to two independent Brownian motions with variance parameter each. Together with 

the appropriate convergence of L^(0), (5.4) then gives the required convergence of H^ , see 
Proposition 5.3 We are now ready to state the main result of this section. 

Theorem 2. For any x > 0, as N ^ 00, 



i7V\ 



^ ^Bl, s > 0}, {Ls{t), s, t > 0}, 



i{Hs,—B^ 

a a 



({Ff , Mi'^, m2'^, s > 0}, {Lf (t), s, t > 0}, S, 

for the topology of locally uniform convergence in s and t. and B^ are two mutually 
independent standard Brownian motions, H solves the SDE (2.1) whose driving Brownian 
motion B is given as 



L is the semi-martingale local time of H, and Sx has been defined in (1.2). 

An immediate consequence of this result is 
Corollary 5.1. For any x > 0, as N ^ 00, 

({/7f ,m1'^,m2'^, s > 0},{L^N{t),t> 0}) {{Hs, ^Bl, ^Bl s > 0},{LsJt),t> 0}) 



a 



a 



in C([0,oo)) X (Z)([0,oo)))= 



Recall (see Lemma 7.1 in the Appendix) that convergence in D([0,oo)) is equivalent to 
locally uniform convergence, provided the limit is continuous. Also note that in the absence 
of reflection, the weak convergence of H^ to Brownian motion would be a consequence of 
Theorem 7.1.4 in |5], and would be a variant of "Rayleigh's random flight model" , see Corollary 
3.3.25 in [20]. 

A first preparation for the proof of Theorem [2] is 

Lemma 5.2. The sequence {H^} is tight in C([0,oo)). 
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Proof: To get rid of the local time term in (5.4), we consider a process R'^ of which H^ 
is the absolute value. More explicitly, let {R^ , W^) be the M x {—1, l}-valued process that 
solves the system (which is exactly (5.3) without reflection) 



R 



N 



N 



2N 
1 + 2 







_ndP:. 



N 



Jo 



We observe that 



{H'',V'') = i\R''\,sgn{R'')W'') 



yN\-,jrN\ 



Clearly tightness of {R^} will imply that of {H^}, since \H. 



rN 



I < \R 



N 



R^\ for all 



s, t > 0. Now we have 



' Na"^ Na^ Na^ 



N 



By Proposition 7.6 in the Appendix, the sequence {R }n>i is tight, and so is {H } 



N>1- 



Proposition 5.3. Fix x > 0. As N ^ 00, 



a a 

in C([0,oo)) X (L)([0,oo)))'' X [0,00), 



where and are two mutually independent standard Brownian motions, and H solves 
the SDE 

H, = -Bs + JLs(O), s > 0, (5.6) 



with Bs := (l/\/2)(-B^ — B'^), and L(0) denoting the local time at level of H. (Note that B 
is again a standard Brownian motion.) 

Proof: The proof is organized as follows. Step 1 establishes the weak convergence of 
(/f^, M^'^, M^'^, L^(0)) along a subsequence. Step 2 and step 3 together characterize the 
law of the limiting two-dimensional martingale, step 4 identifies the limit of the local time 
term. In step 5 we note that the entire sequence converges. Finally step 6 takes the limit in 
the quintuple (including S^). 



Step 1. Note that 
i) from Lemma 5.2 



11) SUPs>0 7^ 



the sequence {H^ , > 1} is tight in C([0, 00]); 
in probability as A — )• 00; 



iii) from Proposition [T^ {M'^'^,N > 1} and {M'^''^,N > 1} are tight in i:>([0,oo]), any 
limiting martingales and being continuous; 



iv) it follows from the first 3 items, (5.4) and Proposition 7.5 that {L^ (0), A > 1} is tight in 
D([0, 00]), the limit K of any converging subsequence being continuous and increasing. 
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Working along a diagonal subsequence we can extract a subsequence, still denoted as an abuse 
like the original sequence, such that along that subsequence 

{H^,M^'^,M^'^,L^iO)) {H,M\M\K). 

Step 2. We claim that for any s > 0, 



in probability, as — t- 00. This follows by taking the limit in the sum and the difference of 
the two following identities : 



rN-_^\dr 



rs rs 

/ l{vN=i}dr+ / l|y 

■70 JO 

Jo Jo 
since H^ /N — )• in probability, as — )• 00, thanks to Lemma 



1}' 

1 uN 



Step 3. By Step 1 iii), Ms = ( ^2 ) > weak limit of ( j^2,Af 1 along the chosen subse- 
quence, is a 2-dimensional continuous martingale. In order to identify it, we first introduce 
some useful notation. We write for the 2x2 matrix whose (z, j)-entry equals MJ x M|, 
and {{M))s for the 2x2 matrix-valued predictable increasing process which is such that 

Mf2-((M)), 

is a martingale, and note that the (z,j)-entry of the matrix {{M))s equals {M^,M^)s. We 
adopt similar notations for the pair 

From Step 2 we deduce that, as — 00, 

-^sl 

in probability, locally uniformly in s, where / denotes the 2x2 identity matrix. Consequently 



'Myy' _^fM}Y" 2 



in I5([0, cxd); M^) as — t- 00, and since the weak limit of martingales is a local martingale, 
there exist two mutually independent standard Brownian motions and such that 



= —Bt, = —Bt, s>0. 
a a 
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Taking the weak limit in (5.4) we deduce that 

h. = ^{bI-b',) + Ik, 

a Z 

= + \k., 

a I 

where B^ = (Bl — i?^)/v2 is also a standard Brownian motion. 



Step 4. For each £ > 1, we define the function fi : M+ — > [0, 1] by fi{x) = (1 — lx)~^ . We 
have that for each A^, £ > 1, s > 0, since L^(0) increases only when H^ = 0, 



Thanks to Lemma 7.3 in the Appendix we can take the limit in this last inequality as N ^ 00, 
yielding 

fe{Hr)dKr - K,^ >0. 

Then taking the limit as ^ — t- 00 yields 

^[j\{Hr=o}dKr-Ks) >0. 

But the random variable under the expectation is clearly nonpositive, hence it is zero a.s., in 
other words 



Ks= l{Hr=0}dKr, Vs > 0, 



which means that the process K increases only when Hr = 0. 
Prom the occupation times formula 



4 



g{Hr)dr= / g{t)L,{t)dt 
Jo 



applied to the function g{h) = l|/i=o}i deduce that the time spent by the process H at 
has a.s. zero Lebesgue measure. Consequently 



hence a.s. 



/ l{Hr=o}dBr = a.s. 
Jo 

Bs= [ l{Hr>0}dBr Vs>0. 

Jo 



It then follows from Tanaka's formula applied to the process H and the function h ^ that 
K = L{<d). 

Step 5. We have proved so far that Q along some subsequence, where = 

{H^, Mi'^, M2'^, L^(0)), Q = {H, ^B^, ^B^, L(0)). Note that not only subsequences but 
the entire sequence Q^,Q^,Q^, . . . converges, since the limit law is uniquely characterized. 
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Step 6. It remains to check that for any x > 0, as — )• 00, 

(Q^, 5f ) {Q, S,) in C([0, 00]) x {Di[0, 00]))' x [0, 00]. 
To this end, let us define the function <I> from M+ x C^(]R+) into M+ by 

4 

y) = inf{s > : y(s) > -i^x}. 

For any x fixed, the function ^{x^ .) is continuous in the neighborhood of a function y which 
is strictly increasing at the time when it first reaches the value x. Clearly Sx = $(x,L.(0)). 
Define 

S'f:=^x,L^m. 

We note that for any x > 0, s 1— t- Ls{0) is a.s. strictly increasing at time Sx, which is a 
stopping time. This fact follows from the strong Markov property, the fact that Hs^ = 0, 
and ie(0) > 0, for all e > 0. Consequently Sx is a.s. a continuous function of the trajectory 
L (0) , then also of Q, and 

(Q^,5;^)^(Q,5.). 

It remains to prove that S'J^ — S^ — )• in probability. For any y < x and N large enough 

— '^x X — X y ' 
Clearly S'J^ — S'J^ =^ Sx — Sy, hence for any e > 0, 

< limsupP(5l^ - > e) < limsupP(5^^ - S'J^ > e) < F{Sx -Sy>e). 

N N 

The result follows, since Sy — )• Sx- as y ^ x,y < x, and Sx- = Sx a.s. ■ 
For the proof of Theorem [2] we will need the following lemmata: 

Lemma 5.4. For any s > 0, t > 0, the following identities hold a.s. 

{H^ - t)+ = 2N r V;^l|^A.>,}dr, 
Jo 

a^N 



Jo 



Proof: The first identity is elementary, and is true along any piecewise linear, continuous 
trajectory {H;^} satisfying dH^ /ds = 2NV^ for almost all s, with G { — !> !}• The other 
identities which we will state in this proof are true a.s. In these identities we exclude the 
trajectories of H^ which have a local maximum or minimum at the level t. This implies that 
the two processes s — t- and s — >• l|j|^jv>j| do not jump at the same time. Hence from 

l{Hiv>t} = l{//,^=t}K^ - l{v;iv=_i}l{j/iv=t} 

0<r<s 

a'^N 



2 

we deduce by differentiating the product that 



fv.^dL^it), 
Jo 



Since (Vj^)"^ = 1, this is the second identity in the lemma. 
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Lemma 5.5. Denote by Ls{t) the local time at level t up to time s of H. Then with probability 
one {s,t) I—)- Ls{t) is continuous from R_|- x into M. 



Proof: This is Theorem VI. 1.7 page 225 of Revuz, Yor [19]. ■ 

Proposition 5.6. For each d > I, < ti < t2 < ■ ■ ■ < t^, 

{ (i^f , Lf (ti), Lf (t2), ..,L^{td)) ,s>0}^{ {Hs, Ls{ti),Ls{t2), .., Ls(td)) , s > O} 

in C([0,oo)) X {D{[0,oo))y. 

Proof: We prove the result in case d = 1 only, the proof of the general case being very 
similar. From (5.3) and Lemma 5.4 we deduce that for any t >0, a. s. 

Lf (t) = 2{H^ - t)+ + J^vfl^H^^t} - 2 r l{H->t}(dMi'^ 

J 

Let 



dM, 



2,N\ 



s > 0. (5.7) 



Jo 



By Proposition 7.6 we have that {U'^}]\f>i is tight in Z)([0,oo)). Moreover, 



M 



2,N\ 



Is — 9 ■ 



Jo 



From the occupation times formula 

'i-{Hr=t}dr = — l{r=t}Ls{r)dr = a.s. 



dr. 



(5.8) 
(5.9) 



Then by Lemma 7.2 from the Appendix we deduce that along an appropriate sequence 



^{HN>t}dr,s> o| 



^{Hr>t}dr,s> o| 



(5.10) 



From (5.8), (5.9) and (5.10), we have again along an appropriate subsequence 

{'^ J\{H^^,,ydBr,'^Bs) in {D{[0,oc)))\ 



N ,^1,7V 



(C/f,M, 



M. 



2,N\ 



Moreover, arguments similar to that used in the proof of Proposition 5.3 establish that 

{H^,U^) {H,,'^ ri^H^^tydBr) in (Z?([0,oo)))l 
^ Jo 

Now from any subsequence, we can extract a subsequence along which we can take the weak 
limit in (5.7). But Tanaka's formula gives us the identity 

Lsit) = 2{Hs -t)+- - j l{H^>t}dBr, 

<^ Jo 

which characterizes the limit of as the local time of H. Since the law of H is uniquely 
characterized, the whole sequence converges. ■ 
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Proposition 5.7. For each s >0 fixed, {L^ > 0}Ar>i is tight in D{[0,oo)). 
Proof: We have 

Lf (t) = 2{Hf - tr + ^y,^l|Hiv>t} + K^l|^,,.>,|dM,^ (5.11 



t I 



where 



K[^ = 2iH^-ty + J^,Vfl,^.,,„ 



From 



{< = 2i/f + — ,y,^l^^^.>ol, > 1} is tight and 



2 

hmsup l/Cf -Klf\< 2\t - t'\ 



it foUows from Theorem 15.1 in [3j that the sequence {/C.^Itvm is tight, and any limit of a 
converging subsequence is a. s. continuous. 



7.4 



We next show that the sequence {G^}Ar>i satisfies the conditions of Proposition 
Condition (1) follows easily from the fact that E(|G^p) < 16s/cj^. In order to verify 
condition (2), we will show that for any T > 0, there exists C > such that for any < t < 

In order to simplify the notations below we let 

tpj. .- l/r---'-{t-e<H^f<t}; 

An essential property, which will be crucial below, is that (p^ ^ijj? = 0- Also (9?^)^ = l^^^l, and 
similarly for if:^ , since those functions take their values in the set { — 1,0, 1}. The quantity 
we want to compute equals up to a fixed multiplicative constant 

'-Jo Jo 
We note that we have the identity 

{f\'^dM?f = 2 r r ip^dM^^'^dM^ + r \^'^\dM^ + a^N^ n^rldr, 
Jo Jo Jo Jo Jo 

and similarly with 93^ replace by ip^ . Because tp^ tp^ = 0, the expectation of the product of 

f r ^^dM^^'^dM^ or r\v'?\dM^ 
Jo Jo Jo 
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with 



JO 



dM^ dMi" or 



N 



\i>f\dM^ 



vanishes. We only need to estimate the expectations 

E 



/o Jo 

and ^ \^r\dr ^ |dr) 







together with similar quantities with ip^ and ip^ interchanged. The estimates of the first 
two expectations are very similar. We estimate the second one as follows, using the Cauchy- 
Schwarz inequality, and Lemma 5.8 below : 

E 



^^\dM^ j \^P^\dr^<^Ej \ip?\d{MN)r^E^(^j \i^^\dry 



Finally, again from Lemma 5.8 



E 



dr l^f |dr) < Ce^ 



The first quantity should be multiplied by A^^, and the second by A^^, and then both should 
be divided by A^*^. The proposition now follows from Proposition 



7.5 



Lemma 5.8. Let s,e,T > 0. Then there exists a constant C such that for all N > 1 and 
0<t,t' <T, 



^{t-e<H?<t}dr) < Ce, 

^{t'~e<H^<t'}dr^ < Ce^ 



^{t-e<HN<t}dr 



Proof: We will prove the second inequality, the first one follows from the second one with 
t = t' and the Cauchy-Schwarz inequality. 

For s,t > define F^{t) := l^Q^j^N^^ydr. It follows readily from the definition of 
that 



~dt 



N 



-it) 
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Hence 



= ^^¥.[f^ J'^ L^{r)L^{u)drdu) 



a 
16 

A 



4 rt ft' 

/ ¥.{L^ {r)L^ {u))drdu 
t-e Jt'-e 



<^e' sup E(Lf(r)LfH) 

■LO 0<r,u<T 

= ^e' sup E((Lf(r))2). 

■LO 0<r<r 



On the other hand, since by Ito's formula there exists a martingale such that 
we conclude that 

supE((Ff )2) < 00. 

Af>l 



The second inequality now follows from (5.11). ■ 

Proposition 5.9. For all d > I, < si < S2 < ■ ■ ■ < Sd, 

(F^,L^,...,Lfj^(^,L,,,...,L,J m C([0,oo))x(Z)([0,oo)))'^. 

Proof: We prove the result in the case d = 1 only, the proof in the general case being very 
similar. From Proposition 5.6 there follows in particular that for all /c > 1,0 < ti < t2 < 
■ ■ • < tfc , we have 

(i/^,Lf(ti),Lf(t2),..,Lf(tfc)) {H,Ls{h),Ls{t2),..,Ls{tk)) 

That is, {L^} converges in finite-dimensional distributions to {Lg}, jointly with iJ^. By 



Proposition 5.7 {L^{t),t > 0}Ar>i is tight. The result follows. ■ 

We are now prepared to complete the 
Proof of Theorem [2j The main task is to combine the assertions of Propositions 5.6 and 



5.9, which means to turn the "partial" convergences asserted for L in these propositions 



into a convergence that is joint in s and t. We will also combine this result with Proposition 



5.3 in order to get joint convergence of all our processes. To facilitate the reading, we will 



divide the proof into several steps. 

Step 1. Let {sn,n > 1} denote a countable dense subset of M+. Our first claim is that for 
ah n G N, 

{H^, M^'^, M^'^, ) {H, Ib^L,,,..., , 5.) (5.12) 



in C(M+) X D{ 



\n+2 
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To make the core of the argument clear, let us write just for the moment 



Y 



N 



(M 



1,N 



a a 



N 



Then (5.12) translates into 



{H,Y,K). 



By Proposition 



5.3 



N vN\ 



(H'\Y 



{H,Y), and by Proposition |5.9| (iJ^,A^^ 



(5.13) 
{H,K). Be- 



cause in our situation A is a.s. a function of H, these two convergences imply (5.13). (More 
generally, this implication would be true if Y and A would be conditionally independent given 
H.) 



Step 2. Now having established (5.12), it follows from a well known theorem due to Skorohod 
that all the processes appearing there can be constructed on a joint probability space, such 
that there exists an event M with P(AA) = and for all uj ^ A/", 



(5.14) 



{H^{u),M}'^{u),M^'''{u)) 
and for all n > 1, 



2 2 

{Hs{u}), —BI{u!), —B'^{uj)) locally uniformly in s > 0, 
a a 

(5.15) 



L^^{t){uj) — Ls„{t){uj) locally uniformly in t > 0, 



(5.16) 



as iV — 7- 00. Here we have made use of Lemma 5.5, which allows us to assume that (s,t) 1— 
Ls{t){ijj) is continuous from M_|_ x M+ into M for all oj ^ AA, possibly at the price of enlarging 
the null set M, and of Lemma 7.1 from the Appendix. 

Step 3. We claim that in the situation described in the previous step one even has for all 

c, r > 0, w ^ AT, 

sup \L^{t,oj) - Ls{t,uj)\^{), (5.17) 

0<s<C,0<t<T 

as N ^ 00. In other words, in Skorokhod's construction there is a.s. convergence of {t) to 
Ls{t), locally uniformly in s and t. To prove ( |5.17| ), we will make use of the fact that for any 
u! ^ Af, and all N,t, the mapping s 1— )• L^{t){io) is increasing and the mapping s 1— )• Ls{t){uj) 
is continuous and increasing. Moreover, since the mapping {s,t) 1— t- Ls{t,uj) is continuous 
from the compact set [0, C] x [0,r] into M+, for any e > 0, there exists 6 > such that 
< s < s' < C, < t < T and s' - s < (5 implies that 

Lsi{t,uj) - Ls{t,u) < e. 

Hence there exists k > 1 and =: sq < ri < • • • < := C such that {rj,0 < i < k] d 
{sn, n>l} and moreover, rj — rj_i < 6 for all 1 < i < A;. We have 



sup \Lg {t,uj) - Ls{t,u)\ < sup [An^i + Bna 

Q<s<C,0<t<T l<i<k 



where 



An,! 



sup {L^{t,io)-Ls{t,uj)y 

ri-i<s<rifl<t<T 

sup (Lf {t, Oj) - Ls{t, U}))' 

ri-i<s<ri,0<t<T 
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For rj_i < s < Tj, 



<{L^^{t,u)-Lr,{t,u:))+ + e, 
iL^{t,io)-Ls{t,u;))- < (L,^_^(t,a;)-L,(t,w))- 

<{L^^_^{t,u)-Lr,_,{t,u;))- +e. 



Finally, 



sup \Lg {t,uj) — Ls{t,uj)\ < 2 sup sup (t, w) — L^. (t, c<j)| + 2e, 

0<s<C,0<t<T 0<i<fcO<i<T 



while from (5.16), 



limsup sup \Lg {t,uj) — Ls{t,uj)\ < 2e. 

AT-s-oo 0<s<Cfl<t<T 



This implies (5.17), since e > is arbitrary. The assertion of Theorem [2] is now immediate 
by combining (5.14), (5.15) and (5.17). ■ 



6 Change of measure and proof of Theorem [T] 

As in the previous section, let, for fixed G N, H^ be a process that follows the dynamics 
described in Proposition |3 . 1 1 for 9 = ^ = 0. We denote the underlying probability measure by 
P, and the filtration by = (J-s)- Our first aim is to construct, by a Girsanov reweighting of 
the restrictions a measure under which H^ follows the dynamics from Proposition 



3.1| for a prescribed ^ > and 7 > 0. 

Here, a crucial role is played by the point process of the successive local maxima 
and minima of H'^ , excluding the minima at height 0. Under P, this is a Poisson process 
with intensity a'^N'^. More precisely, the process Q^''^ which counts the successive local 
minima of H^ (except those at height 0) is a point process with predictable intensity A^'^ — 

4.1 j;„j.„l,l„ ;„i„„„;j.„ \ 2,Ar Ar2_2-i 

^{V;^=+l}- 

is the (cadlag) sign of the slope of H'^ .) 

For the rest of this section we fix > and 7 > 0. In view of Proposition |3 . 1 1 we want to 



A^^(T^l|yjv^_]^|, and the process Q^'^ which counts the successive local maxima of H^ is a 
point process with predictable intensity Ag'^ := -/V^cr^lryiv^ , ^1 . (Recall that the process V'^ 



change the rate Xl'^ to As'^(l + ■^^) and the rate As'^ to Xl'^ {1 + '^^^J^^"-* ). 

As in Section [5] we will use the process = P^ — Na'^s, s > 0, which is a martingale 



under P. Taking the route designed by Proposition |7.13| in the Appendix, we consider the 
local martingales 

" '~ Jo iV^ {^r^^ ' ' N ^{v;^=i}^^^ ' ^ .-X '+X '. 

Let Y'^ := £{X^) denote the Doleans exponential of X^ . Proposition 7.9 in the Appendix 
recalls this concept and the fact that Y'^ is the solution of 

Yf = l + J\jl^-^^l,y.^_,, + :^^^^ s>0. (6.1) 



We will show that Y^ is a martingale, which from Proposition 7.13 will directly render the 
required change of measure. 
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Proposition 6.1. Y is a {T ,¥) -martingale. 

Proof: Under P, is a positive super-martingale and a local martingale. It is a martingale 
if and only if 

E[y/] = yj^ = 1, (6.2) 

which we will show. A key idea is to work along the excursions of H^ , that is, along the 
sequence of stopping times ri^'* := S^i^ As, a = 0, 1, 2, . . . Since and s are fixed, we will 

suppress the superscripts N and s for brevity and write Tq instead of Ta'^ ■ 



Step 1. We want to show that 



E[Y^n = l, a = 1,2,... (6.3) 



Between Ta~i and r^, the solution of (6.1) is bounded above by the solution of the same 



equation with replaced by , which takes the form dYj^ = Y^a^ dP^ . If we denote 



by {Tfc, A; > 1} the successive jump times of the Poisson process P , we have or each k 
Y^ = Yr^_{l + tty ). Consequently for all Ta-i < t < Ta, 

Yr"" < n (1+ 2' 1 I ^^tSK\ ^ 

^<^-l k>l:Ta-l<Tk<Ta 

Within the excursion of H^ between the times Ta-i and Ta, these jump times coincide with 
the times of the local maxima and minima of H^ in the time interval {Ta-i,Ta). Since for 
a > 1 there are reflections of at in the time interval (0,ra_i), the parity of those k for 
which V^^_ = 1, Ta-i <Tk<Ta, depends on a. However, noting that 

^T,(^£)<^A;, (6.4) 
we infer the existence of a constant c > such that 

yN 

- <C^- X(P,^ + 1)!!, Ta^l<r<Ta, (6.5) 



where we define for k G N, A;!! = 1 • 3 • 5 • • • A: if A; is odd, and A!! = 1 • 3 • 5 • • • (A — 1) if A is 
even. 
Now 

/ f-Ta yN \ 

= y^l^ 1 + / ^[p(r) + L^_{H^)q{r)]dM^'^ , (6.6) 

where < p{r) < §A^, < q{r) < jj, p and q are predictable. The claimed equalities 
E[y^^] = l,o = l,2,..., fohow by induction on a from (|6.6l), provided the process 

rs yN 

Mf= l].„_,.,](r)-^b(r)+L^_(i?f)g(r)]dM,^, . > 0, 

•''0 

is a martingale. From Theorem T8 in Bremaud [4j page 27, this is a consequence of the fact 
that 



E 



rs yN 

/ hra-ura]ir)TP^\p{r) + Lt[H:')q{r)]dr < oo. 

•''0 
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In order to verify the latter inequality, we compute 



E 



[p{r) + L^_{H^)q{r)]dr 



yN 

Ta-l Ta-1 



< CnsE 



< CnsE 



x(P,^ + l)!!(l + P-) 



c^=" X (Pf + 1)!!(1 + P, 



< CnsCn,s, 



where we have used (6.5), ( |6.4[ ) and Ta < s, and where Cat and Cn^s are constants which 
depend only on N and {N,s), respectively. The fact that Cn^s < oo follows from 



E[c^'='^(Pf + l)!!Pf ] = exp{-A^Vs}^c^(A; + 1)!!A; 



k=0 



Since 



{k + l)l\k _ k{k + l) _kik + l) 1 
2.4...(2[§]) " 2[il 'm 
1 



kl 



< 



[|]! 



Vfc > 20, 



we deduce that E[c2'' {P^ + l)!!Pj^] < oo. This completes the proof of (6.3). 

Step 2. We can now define a consistent family of probability measures on J>^, a 

1,2, ... by putting 



^N,s,a 



Yf^, a e N. 



We write P^'* for the probability measure on the cr-field generated by union of the u-fields 
J^^^ , a = 1, 2, . . ., whose restriction to Jv^ is ]P^>*'" for all a = 1, 2, . . ., and put 

A := inf{a G N : Tq = s}. 

We will now show that 

1. (i) A < oo P^'*- a.s. (and consequently ta = s P^'''- a.s.), 

2. (ii) under P^''^, (-ff/^)o<r<TA = {H^)o<r<s is a stochastic process following the dynamics 



specified in Proposition 3.1 



Indeed, applying Girsanov's theorem (Proposition 7.13 in the Appendix) to the 2-variate 
point process 



1,N r^2,N\ 



{Qr,Q 



we have that under p-'^'*''^ 



Ql^^ has intensity (iVV^ + 26'iV)l 



dP, 



N 



ndP^ ], < r < r„ 



(6.7) 



N/ttN 



dr 



Qp'^ has intensity a'[N^ + jNU; {W; )]l^yN_^^^dr. 
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Thus, for all a e N, (H^) 

o<r<Ta is, under p^'^'"-^ a stochastic process following the dynamics 



{Hr)Ta-l<r<Ta, ^ 

(i) and (ii). 



from Proposition |3.1| up to the stopping time . Considerin g th e sequence of excursions 



1,2,... under 



'■s>N,s 



we infer from Lemma 



3.3 



the validity of the claims 



Step 3. We now prove (6.2). For this we observe that 



E[y/] = Y,ny!';A = a] = Y,nYr1:A = a] = 5^1P^'^(A = a)= P^'^(^ < oo) = 1. 



a>l 



a>l 



a>l 



Corollary 6.2. Let be the probability measure on T whose restriction to J-g, s > 0, has 
density YJ^ (given by (6.1)) with respect to P|j-j. Then under F'^ the process H^ follows the 



dynamics from Proposition\3. 1\ for the prescribed 6 and 7. 



Proof: This is immediate from Proposition 6.1 and the discussion preceding it, combined 



with Proposition 7.13 in the Appendix applied to the process defined in (6.7), now with 
< r < 00. ■ 

Next we will analyze the behaviour of the Girsanov densities as — t- 00. For this we use 
the two martingales M^'^ and M^'^ defined in (5.5), and note that (6.1) can be rewritten as 



TV 



1 + 



JO 2 



The two (pure jump) martingales M^'^ and M'^'^ have jump sizes 2/{Na'^), hence the random 



2,N 



7.10 



variable under the expectation in formula (7.6) vanishes for suitably large N. Thus (see 
in the Appendix), the sequences {M^'^}7v>i and {M^'''^}7v>i have uniformly 



Definition 

controlled variations, and because of Proposition 7.11 (1) they are "good". Hence 



Jo 2 

Moreover, by Proposition 



7.11 



(3), {Xf}i\i>i is also a good sequence, hence by Proposition 



7Tn(2) 



Y 



N 



£{X 



£{X) =: Y. 



Combining these facts with Corollary 5.1, we deduce, again from Proposition 7.11 (3), that 



{H,Ls^,Y). 



(6.8) 



Since and are mutually orthogonal, by Proposition 7.9 we have 
Y,=£[ 



cr 2 



a 



\/2jaLriHr) ^^2 



expj— + 



V2-faLr{Hr) 



dBl- [^^ + ^-^Lr{Hrf]dr] 
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Applying Lemma 2.2 and Lemma 2.3 again, we deduce that y is a martingale. In particular 
E[Ys] = 1 for all s > 0. Define the probability measure P by 



Ys, Vs > 0, 



then H, under P, solves the SDE (1.3) with Bs there replaced by 

1 



B, := ^(Bl - Bl) - -s- 



70" 



Lr{Hr)dr, 



which is a standard Brownian motion under P due to Proposition 7.12 

The following general and elementary Lemma will allow us to conclude the required con- 
vergence under the transformed measures. 

Lemma 6.3. Let {(,n,i1n), iCyV) random pairs defined on a probability space {fl,T,F), 
with r]js!, rj nonnegative scalar random variables, and (,i\f, ^ taking values in some complete 
separable metric space X. Assume that E[r/7v] = IE[r/] = 1. Write {£,n,Vn) for the random 
pair (Cn^Vn) defined under the probability measure which has density rjN with respect to 
P, and {rj, ^) for the random pair (rj, ^) defined under the probability measure P which has 
density rj with respect to P. Then {S,n^Vn) converges in distribution to {fj,^,), provided that 
{Cn,Vn) converges in distribution to {(,,r]). 

Proof: Due to the equality E[r/jv] = ^[n] = 1 s-^d a variant of Scheffe's theorem (see Thm. 
16.12 in [2j), the sequence t/tv is uniformly integrable. Hence for all bounded continuous 
F -.X xR,-^R, 



Combining (6.8) with Lemma 6.3 yields the 



Theorem 3. Let H^ be a stochastic process following the dynamics specified in Proposition 
3.1. and let H be the unique weak solution of the SDE (1.3). We have 



(ttN jN n 



iH,LsJ m C([0,oo]) xZ)([0,oo]), 



(6.9) 



where S^ and Sx are defined in (5.2) and (1.2). 



We can now proceed with the 
Completion of the proof of Theorem 



2;N,x fQ^iQ^g ^}jg dynamics (3.1). From (6.9), ^Ls^ is the limit in distribution of Z 
N ^ oo. Hence by Proposition 



EM ij : 



Define Z, 



'^•":=^i^iv(t). By Corollary 



3.2 



N,x 



as 



4.3 



completes the proof of Theorem [IJ 



1 1—)- ^Ls^{t) is a weak solution of the SDE (1.1), which 
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Remark 6.4. Theorem [T] establishes a correspondence between the solution H of the SDE 



(1.3) and the logistic Feller process, i.e. the solution of (1.1). This connection can be expressed 
in particular through the occupation times formula for H, which states that for any Borel 
measurable and positive valued function /, 

fSx 



/ f{Hs)ds= / f{t)Z^dt. 
Jo Jo 

This formula in the particular case / = 1 states that 

POO 

S,= / Zfdt. 
Jo 

The quantity on the right is the area under the trajectory Z^ . It is the limit of the properly 
scaled total branch length of the approximating forests defined in Section [sj We now 
establish another identity concerning this same quantity, with the help of a time change 
introduced by Lambert in |10j . Consider the additive functional 

At = f Zfdr, 

and the associated time change 







at = inf{r > 0, Ar > t}. 



As noted in |1U| . the process := Z^^ is an Ornstein-Uhlenbeck process, solution of the 
SDE 

dUf = {9- -/Uf)dt + adBt, = x. 

Of course this identification is valid only for < t < r^, where := inf{t > 0, Ut = 0}. Let 
Tx be the extinction time of the logistic Feller process Zf . We clearly have ar^ = T^, and 
consequently 

poo 

Zfdr. 





We have identified the time at which the local time at of the exploration process H reaches x 
with the area under the logistic Feller trajectory starting from x, and with the time taken by 
the Ornstein-Uhlenbeck process to reach 0. The reader may notice that in the particular 
case 7 = 0, the identity Sx = is not a surprise, see also the discussion and the references 
in [15] Section 6. 



7 Appendix 

7.1 Skorohod's topology and tightness in D([0, oo)) 

We denote by -D([0, oo)) the space of functions from [0, oo) into M which are right continuous 
and have left limits at any t > (as usual such a function is called cadlag). We briefly 
write D for the space of adapted, cadlag stochastic processes. We shall always equip the 
space -D([0, oo)) with the Skorohod topology, for the definition of which we refer the reader 
to Billingsley |3] or Joffe, Metivier [7j. The next Lemma follows from considerations which 
can be found in [3], bottom of page 124. 
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Lemma 7.1. Suppose {x„,n > 1} C I?([0,oo)) and Xn ^ x in the Skorohod topology. 

(i) If X is continuous, then Xn converges to x locally uniformly. 

(a) If each Xn is continuous, then so is x, and x„ converges to x locally uniformly. 
In particular, the space (7([0, oo)) is closed in L'([0,oo)) equipped with the Skorohod topology. 



The following two lemmata are used in the proofs of Propositions 5.6 and 5.3 
Lemma 7.2. Fix t > 0. Let Xn,x G C([0,oo)),n> 1 he such that 

1. Xn ^ X locally uniformly, as n ^ oo. 

2. for each s > 0, 



f 

Jo 



^{x{r)=t}dr = 0. 



Then 

/ ^Mr)>t}dr / l{^.(,,)>t}dr locally uniformly in s>0. 
Jo Jo 

Proof: We prove convergence for each s > 0. The local uniformity is then easy. Given 
£ > 0, there exists A'o such that 



Then for all n > Nq, 



sup |x„,(r) — x(r)| < e Vn > A^o- 

0<r<s 



|l{a;„(r)>t} " l{x-{r)>t}l ^ l{i-£<x{r)<t+e} 



PS PS PS 

/ ^{x„lir)>t}dr - / l{^(r)>t}(ir < / l{t-e<xir)<t+e}dr. 

Jo Jo Jo 



The result follows from 



PS PS 

/ ^{t~e<x(r)<t+e}dr = / l{^(,.)=t}dr = 0. 
^~^"J0 Jo 



Lemma 7.3. Let Xn,yn £ ^([0, oo)), n > 1 and x, y € C([0, oo)) be such that 

1. for all n> 1, the function t — )■ yn{t) is increasing; 

2. Xn ^ X and yn — y, both locally uniformly. 
Then y is increasing and 

/ Xn{s)dyn{s) — )• / x{s)dy{s), locally uniformly in t > 0. 
Jo Jo 
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Proof: We prove convergence for each t > 0. The local uniformity is then easy. 



x{s)dy{s) 



Xn{s)dy„ 



< 



[x{s) - Xnis)]dynis] 







+ 



x{s)[dy{s) - dynis)] 



- Us)\[dyis) - dynis)] + / |C,(s)| - 



< sup \x{s) - Xn{s)\yn{t) + 
0<s<t 

where is a step function which is such that supQ<g<^ — £,eis)\ < £• The first and last 
term of the above right hand side clearly tend to as n — )• oo. Then 



lim sup / x{s)dy{s) 



/ Xn{s)dyn{s) < elimsup[y„(t) + 

JO n— >oo 

< 2y{t) X e. 



It remains to let e — >• 0. ■ 

We first state a tightness criterion, which is Theorem 13.5 from [3] : 

Proposition 7.4. Let {X^^ t > 0}n>i be a sequence of random elements of D{[0,oo)). A 
sufficient condition for {X^} to be tight is that the two conditions (i) and (ii) be satisfied : 

(i) For each t > 0, the sequence of random variables {X", n > 1} is tight in M; 

(ii) for each T > 0, there exists (3,C > and 9 > 1 such that 



E X 



H+h 



xiflxi 



X 



t-h\ 



for alio <t <T, < h <t, n> 1. 



Note that convergence in -D([0, oo)) is not additive : x„ — )• x and yn ^ y in -^([0, oo)) does 
not imply that Xn + yn — ^ x + y in D{[0, oo)). This is due to the fact that to the sequence x„ is 
attached a sequence of time changes, and to the sequence yn is attached another sequence of 
time changes, such that the time changed x„ and yn converge uniformly. But there may not 
exist a sequence of time changes which makes Xn + yn converge. If now {X]}, t > 0}ri>i and 
{y^", t > 0}„>i are two tight sequences of random elements of D([0, oo)), we cannot conclude 
that {Xf + y^", t > 0}„>i is tight. However, if and yn ^ y in D{[0,oo)) and x is 

continuous, then we deduce easily from Lemma 7.1 that Xn + yn — ^ x + y in Z)([0,oo)). It 
follows 

Proposition 7.5. If {X^, t > 0}n>i and {Y/^, t > 0}n>i are two tight sequences of random 
elements o/i5([0, oo)) such that any limit of a weakly converging subsequence of the sequence 
{X", t > 0}„>i is a. s. continuous, then {X" + 1^", t > 0}„>i is tight in D{[0,oo)) . 

Consider a sequence {X",t > 0}n>i of one-dimensional semi-martingales, which is such 
that for each n > 1, 

ft 

^niX^)ds + M^, t>0; 



Xf 



t > 0; 
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where for each n > 1,M" is a locally square- integrable martingale, (pn and ipn are Borel 
measurable functions from M into M and M+ respectively. 

The following result is an easy consequence of Theorems 16.10 and 13.4 from [3]. 

Proposition 7.6. A sufficient condition for the above sequence > 0}n>i of semi- 

martingales to be tight in L'([0,oo)) is that both 

the sequence of r.v.'s {X^,n > 1} is tight; (7-1) 

and for some p > I, 

VT > 0, the sequence of r.v. 's { [ [| (pn{X^) \ +'4>n{X]^)fdt,n > l] is tight. (7.2) 

Jo 

Those conditions imply that both the bounded variation parts {V^,n > 1} and the martingale 
parts {M'^,n > 1} are tight, and that the limit of any converging subsequence of {V""} is a.s. 
continuous. 

If moreover, for any T > 0, as n ^ oo, 

supo<t<T I -^i" ~ probability, 

then any limit X of a converging subsequence of the original sequence {X"}„>i is a.s. con- 
tinuous. 



Remark 7.7. A sufficient condition for (7.2) is that for all T > 



{ sup [\MXt)\ + MXt)], n > 1} is tight. (7.3) 

0<t<T 



Remark 7.8. A sufficient condition for (7.2) is that for all T > 



hmsup sup E[^l{X^) + V'(^r)] < oo. (7.4) 

n>l 0<t<T 



Indeed, (|7.4|) yields 

rT 



limsupE / [\ipn{X^)\+MXt)?dt<oo, 
n Jo 



which in turn implies (7.2). 



7.2 Doleans exponential and "goodness" 

For a cadlag semi-martingale X = {Xt,t > 0), consider the stochastic linear equation of 
Doleans ^ 

Yt = l+ [ Yr-dXr. (7.5) 







The following proposition follows from Theorem 1 and Theorem 2 in [12], page 122. 



Proposition 7.9. (1) Equation (7.5) has a unique solution (up to indistinguishability) within 
the class of semi-martingales. This solution is denoted by £{X) and is called the Doleans 
exponential of X . It has the following representation 

£{X)t = exp {Xt -Xo- ^mt} 11(1 + AX,)e-^^^ 

r<t 
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(2) If U and X are two semi-martingales, then 

£{U)t£{X)t = £{U + X + [U,X])t. 

(3) If X is a local martingale, then £{X) is a nonnegative local martingale and a super- 
martingale. 



For (5 > we define hs : M+ — > M+ by hs{r) = (1 — 5/r)^ . For x G L)([0,oo)), we define 
e ^([O,^)) by 

xf := xt - ^ hsd Axs \)Axs. 

0<s<t 

Definition 7.10. (1) Let G,G^ in D,{G",n > 1} be a sequence of semimartingales adapted 
to a given filtration (J^t) o-nd assume G" =^ G as n —t- oo. The sequence (G") is called good 
if for any sequence {I'^,n > 1} of [Tf) -progressively measurable processes in D such that 
(I'^,G"') =^ (IjG) as n ^ oo, then G is a semi-martingale for a filtration with respect to 
which I is adapted, and (/", G", / I^_dG'^) =^ (/, G, / Is-dGg) as n ^ oo. 
(2) A sequence of semi-martingales {G"'}„>i is said to have uniformly controlled variations 
if there exists 5 > 0, and for each a > 0,n > 1, there exists a semi-martingale decomposition 
Qn,5 ^ ]^n,5 _^ j^n,s ^ stopping time T"'" such that P({r"'" < a}) < ^ and furthermore 

supEjiM^'^M^'^jtAT":" + / I d^"''' I I < oo. (7.6) 

n ^ Jo J 

It follows from pages 32 ff. in [9] 

Proposition 7.11. Let G,G^ in D,{G",n > 1} be a sequence of semi-martingales and 
assume G" =^ G. 

(1) The sequence {G"} is good if and only if it has uniformly controlled variations. 

(2) //{G"} is good, then (G",£:(G'^)) {G,£{G)). 

(3) Suppose (I",G") {I,G), and {G"} is good. Then J" = fI2_dG'^, n = 1,2,..., 

is also a good sequence of semi-martingales. Moreover under the same conditions, 
(/",G",f(G")) {I,G,£{G)). 

7.3 Two Girsanov theorems 

We state two multivariate versions of the Girsanov theorem, one for the Brownian and one 
for the point process case. The second one combines Theorems T2 and T3 from [3], pages 
165-166. 

Proposition 7.12. Let {{B^\ . . . ,B^f^), s > 0} be a d-dimensional standard Brownian 
motion defined on the filtered probability space {VL,F,¥). Moreover, let cp = (0i, ...,(j)d) be an 
T -progressively measurable process with (l)i{r)'^dr < oo for all 1 < i < d and s > 0. Let 
X^s^ := /q (l)i{r) dB^'^ and put Y := £{XW -\ ^ X^), or in other words 

Ys = exp { j\cp{r),dBr) -\ \ ? dr]. 
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//E[ys] = 1, s > 0, then := — (f){r)dr, s > 0, is a d- dimensional standard Brownian 
motion under the probability measure P defined by dF / dF = s > 0. 

Proposition 7.13. Let {{Q^s^\ ...,Qi'^^),s > 0} be a d-variate point process adapted to some 
filtration T , and let {Ai*\ s > 0} 6e the predictable (P, -intensity of Q^^\ 1 < i < d. Assume 
that none of the Q^^\ Q^-'\ i 7^ 3, jump simultaneously. Let {iJL^\r > 0},1 < i < d, be 
nonnegative J- -predictable processes such that for all s > and all 1 < i < d 



s 

fi^fhi'Ur < oo P -a.s. 

For i = 1, . . . , d and s > define 

X« := - l)dMW, y« := y = + . . . + X^^^). 

Jo 

Then, with {T^, k = 1,2 . . .} denoting the jump times of Q'^^\ 



y(i) 



n /^S^]) exp { j\l - 4^))A«dr} and y, = J] ^i'^ ^ > 0- 



//" E[ys] = 1, s > 0, i/ien, /or each 1 < i < d, the process Q*-*-* has the -intensity 
X^^ = r > 0, where the probability measure P is defined by dF \jr^/dF = Yg, s > 0. 

Acknowledgements. We thank J.F. Delmas and Ed Perkins for stimulating discussions, 
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